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Abstract 

For applications regarding transition prediction, wing design and control 
of boundary layers, the fundamental understanding of disturbance growth in 
the flat-plate boundary layer is an important issue. In the present work we 
investigate the stability of boundary layer in Poiseuille flow. We normalize 
pressure and time by inertial and viscous effects. The disturbances are taken 
to be periodic in the spanwise direction and time. We present a set of linear 
governing equations for the parabolic evolution of wavelike disturbances. 
Then, we derive modified Orr-Sommerfeld equations that can be applied in 
the layer. Contrary to what one might think, we find that Squire's theorem 
is not applicable for the boundaiy layer. We find also that normalization by 
inertial or viscous effects leads to the same order of stability or instability. 
For the 2D disturbances flow (6 = 0), we found the same critical Reynolds 
number for our two normalizations. This value coincides with the one we 
know for neutral stability of the known Orr-Sommerfeld equation. We no- 
ticed also that for all overs values of k in the case 9 = correspond the same 
values of Res at a = whatever the normalization. We therefore conclude 
that in the boundary layer with a 2D-disturbance, we have the same neu- 
tral stability curve whatever the normalization. We find also that for a flow 
with hight hydrodynamic Reynolds number, the neu- tral disturbances in the 
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boundary layer are two-dimensional. At last, we find that transition from 
stability to instability or the opposite can occur according to the Reynolds 
number and the wave number. 

1 Introduction 

Boundary-layer theory is crucial in understanding why certain phenomena occur. 
It is well known that the instability of boundary layers is sensitive to the mean 
velocity profile, so that a small distortion to the basic flow may have a detrimental 
effect on its stability. Prandtl (1904)[9j] proposed that viscous effects would be 
confined to thin layers adjacent to boundaries in the case of the motion of fluids 
with very little viscosity i.e. in the case of flows for which the characteristic 
Reynolds number, R e , is large. In a more general sense we will use boundary-layer 
theory (BLT) to refer to any large-Reynolds-number. Ho and Denn studied low 
Reynolds number stability for plane Poiseuille flow by using a numerical scheme 
based on the shooting method. They found that at low Reynolds numbers no 
instabilities occur, but the numerical method led to artificial instabilities. Lee and 
Finlayson used a similar numerical method to study both Poiseuille and Couette 
flow, and confirmed the absence of instabilities at low Reynolds number. R. N. 
RAY et al IflOll investigated the linear stability of plane Poiseuille flow at small 
Reynolds number of a conducting Oldroyd fluid in the presence of magnetic field. 
They found that viscoelastic parameters have destabilizing effect and magnetic 
field has a stabilizing effect in the field of flow but no instabilities are found. 

In this paper, we study the linear stability of boundary layer in a plane 
Poiseuille flow. For this, we derive two fourth-order equations name modified 
fourth-order Orr-Sommerfeld equations governing the stability analysis in bound- 
ary layer for the flow. The first is obtained by making dimensionless quantities 
by the inertial effects. The second takes into account the form adopted by the 
rheologists i.e. make the quantities dimensionless by normalizing by the viscous 
effects. This allowed us to see the effect of each type of normalization on the 
stability in the boundary layer. So, we solve numerically the corresponding eigen- 
values problems. We employ Matlab in all our numerical computations to find 
eigenvalues. 

The paper is organized as follows. In the second section the boundary layer 
theory is presented. In the third section we present the general formulation, high- 
lighting the fundamental equations that model the flat-plate boundary layer flow 
according to the normalization by inertial and viscous effects. In the fourth sec- 
tion the modified Orr-Sommerfeld equations governing the stability analysis in 
boundary layer are checked and in the fifth section, analysis of the stability is 
investigated. The conclusions and perspectives are presented in the final section. 
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2 Boundary Layer Theory 



When applying the theory of complex potential around an airfoil in consider- 
ing the model of inviscid incompressible irrotational plan flow, we know that the 
model allows to deduce the lift but the drag is zero. This is contrary to experimen- 
tal observations which show that the drag affects all flows of real fluids. These are 
viscous. They adhere to the walls and the tangential component of the velocity is 
zero if the wall is fixed. The latter condition can be satisfied by the perfect fluid. 
Moreover, the irrotational condition is far from reality as we know that the pro- 
duction of vorticity occurs at the walls. To remedy the deficiencies of the theory 
of perfect fluid, it must appeal to the theory of the boundary layer which is a nec- 
essary correction for flows with high Reynolds numbers. Theoris the boundary 
layer is due to L. Prandtl|9]|. The boundary layer is the area of the flow which 
is close to the wall or of an obstacle present in a uniform flow at the upstream 
infinity or on the confining walls of internal flow. Within the boundary layer is 
a thin zone, it is estimated that viscous effects are of the same magnitude as the 
inertial effects. The boundary layer is the place of intense generation of vorticity 
which will not diffuse into the area outside thereof. This leads to a very modern 
concept of comprehensive approach to the problem by breaking it down into two 
areas: firstly the boundary layer where we will consider the viscous effects in a 
model of simplified Navier-Stokes and other from the outer area where we will 
use the complex potential theory in the inviscid incompressible flow. This outer 
zone has speeds which are of the same order of magnitude as that of the incident 
flow. 

The boundary layer along an obstacle is therefore thin since the fluid travel 
great distances downstream of the leading edge during the time interval during 
which the vortex diffuse only a small distance from the wall. The creation of 
vorticity in the boundary layer allows the physical realization of the fluid flow 
around the profile. This movement gives rise to a wake in the area near the trailing 
edge. The importance of the wake depend on the shape of the obstacle and the 
angle of incidence of the upstream flow at the leading edge. 

We consider incompressible flow of a fluid with constant density p and dy- 
namic viscosity fi, past a body with typical length L. We assume that a typical 
velocity scale is U, and the Reynolds number is given by 

R c = ^ » 1. (1) 

For simplicity we will, for the most part, consider two-dimensional incom- 
pressible flows, although many of our statements can be generalised to three- 
dimensional flowsand/or compressible flows. 

Boundary Layer Theory applies to flows where there are extensive inviscid 
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regions separated by thin shear layers, say, of typical width 5 <C L. For one such 
shear layer take local dimensional Cartesian coordinates x and y along and across 
the shear layer respectively. Denote the corresponding velocity components by u 
and v respectively, pressure by p and time by t. On the basis of scaling arguments 
it then follows that 



i 



(5~i? e 2 L<L. (2) 

Further, it can also be deduced that the key approximations in classical Bound- 
ary Layer Theory are that the pressure is constant across the shear layer i.e. 

= ~Py, (3) 

and that streamwise diffusion is negligible, i.e. if • represents any variable 

*yy ^ > *xx (4) 

The former approximation is more significant dynamically. 
Now, using the transformations 

(x, y, i, u, v,p) -»■ (Lx, R e 2 Ly, U^Lt, Uu,R e 2 Uv, pU 2 p) (5) 

and taking the limit R e — > oo, the Boundary Layer Theory equations can be 
deduced from the Navier-Stokes equations: 

U t + UU X + VUy = + Uyy, (6) 

= ~Py, U X + Vy = (7) 

For flow past a rigid body the appropriate boundary conditions are 

u = v = on y = and u — > U(x, t) as y — > oo, (8) 

where U (x, t) is the inviscid slip velocity past the body. Further, from (6) 
evaluated at the edge of the boundary layer 

-p x = U t + UU x . (9) 
We define the viscous blowing velocity out of the boundary layer to be 

v b (x, t) = lim (v + U x (x, t)y). (10) 

y— >oo 

Vb indicates the strength of blowing, or suction, at the edge of the boundary 
layer induced by viscous effects. It is good diagnostic for dynamically significant 
effects within the boundary layer much better than, say, the wall shear u y (x, 0, t) 
which can remain regular while Vb(x, t) becomes unbounded. 
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3 General formulation 



Consider an incompressible boundary layer over a flat plate as illustrated in the 
following Figure||8]]. 




Pi g u r e : Fl a1>pl ate b on n rl ary 1 aye r flo w . 



The streamwise coordinate x is scaled with the length scale I, which is a fixed 
distance from the leading edge. The wall-normal and spanwise coordinates y 
and z, respectively, are scaled with the boundary-layer parameter 5 = a/ vT/U^ , 
where v is the kinematic viscosity of the fluid and is the streamwise freestream 
velocity at the distance / from the leading edge. The streamwise velocity u 
is scaled with , while the wall- normal and spanwise velocities v and w, 
respectively, are scaled with U^S/l. The pressure p is scaled with pU^S 2 /I 2 
and the time t is scaled with l/U^ . The Reynolds numbers used here are de- 
fined as Rei = Uool/v and Reg = U^S/v. It is useful to note the relations 
1/8 — Res = y/Rei . We want to study the linear stability of a high Reynolds 
number flow. The non-dimensional Navier-Stokes equations for an incompress- 
ible flow where pressure and time are normalized by inertial effects 

du dv dw 

7r + 7r + 7T = ' < n ) 

ox ay oz 

Du 1 dp ^ 1 d 2 u ^ d 2 u ^ d 2 u 
Dt Re 2 dx Re 2 dx 2 dy 2 dz 2 ' 

Dv dp 1 d 2 v d 2 v d 2 v 

— = — - H 1 1 , (13) 

Dt dy Re 2 dx 2 dy 2 dz 2 ' 
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Dw dp 1 d 2 w d 2 w d 2 w 

= — — H 1 1 , (14) 

Dt dz Re 2 dx 2 dy 2 dz 2 ' 

where # = l + ^+^ + ^f 

are linearized around a two-dimensional, steady base flow (U (x, y), V(x, y), 0) 
to obtain the stability equations for the spatial evolution of three-dimensional, 
time-dependent disturbances (u(x, y, z, t), v (x, y, z, t),w(x, y, z, t),p(x, y, z, t)) . 
The base flow and the disturbances are scaled in the same way. The disturbances 
are taken to be periodic in the spanwise direction and time, which allows us to 
assume solutions of the form 

PX 

f = f(x,y)exp(iRe$ / a(x)dx + i/3z — iut) (15) 

Jxo 

where / represents either one of the disturbances u, v, w or p. The complex 
streamwise wave number a captures the fast wavelike variation of the modes and 
is therefore scaled with 1/5. a itself is assumed to vary slowly with x. Since x 
is scaled with /, the factor Res appears in front of the integral. The x-dependence 
in the amplitude function / includes the weak variation of the disturbances. The 
real spanwise wave number (3 and angular frequency to are scaled in a consistent 
way with z and t, respectively. Introducing the assump- tion (1 1) in the linearized 
Navier-Stokes equations and neglecting all third order terms in 1 / Res or higher, 
we arrive at the parabolized stability equations in boundary-layer scalings 

u x + iResau + v y + if3w = 0, (16) 



Px 1 Otp 9 

(U x + iResoiU — iu)u + Uu x + Vu y + U y v + — ^ + — — = u yy — k u, (17) 

Re$ Res 



(V y + iResaU — iuj)v + Uv x + V x u y + Vv y + p y = v yy — k 2 v, (18) 

(iResaU — ico)w + Uw x + Vw y + i(3p = w yy — k 2 w, (19) 
where k 2 = a 2 + (3 2 . 

If we normalize pressure and time by viscous effects i.e. t scaled with pi 2 / fi 
and p scaled with lj the Navier-Stokes equations take the following forms 

du dv dw 

IT + 7T + IT = °> ( 2 °) 
ox ay oz 
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1 du „du A 9« „9m 1 dp 1 d u d u d u 



' +V— +W— = -TT9— + TT97T^ + TT^ + TT"^ ( 21 ) 



i?e^ dt dx dy dz Re 2 dx Re 2 dx 2 dy 2 dz 2 



1 dv ^dv ^dv „dv dp 1 d 2 v d 2 v d 2 v 

+ u— + v— + w— = -^ + —^— + — + —, (22) 



Re 2 dt dx dy dz dy Re 2 dx 2 dy 2 dz 2 



1 dw A dw „ dw „ dw dp 1 d 2 w d 2 w d 2 w 

+ u— + v— + w— = - 7 f + —^— T + — + ^—. (23) 



Re 2 dt dx dy dz dz Re 2 dx 2 dy 2 dz 2 

The linearized Navier-Stokes equations with the previous disturbances under 
the same considerations become 

u x + iResait + v y + if3w = 0, (24) 



iuj p x tap 9 

(U x + iRe s aU - -^-n)u + Uu x + Vu y + U y v + —j + — — = u yy - k u, (25) 
R&fi R^s R&s 



(V y + iRe s aU - -=-~)v + U v x + V x u y + Vv y + p y = v yy - k 2 v, (26) 
Re s 



(iResaU — FTij)^ + Uw x + Vw y + i/3p = w yy — k w, (27) 
Re s 

where k 2 = a 2 + f3 2 . 



4 Modified Orr-Sommerfeld Equation 

Considering temporal the stability problem with u = ac and (a, (3) real, we will 
simplify the problem. Our strategy will be first to eliminate u, w, p to leave a 
single equation in v. This can be used finally to determine the linear stability (or 
instability) in the boundary layer oy the base flow. Remember that c = c r + ici 
and if q < 0, the flow is stable; q > 0, the flow is unstable and we have neutral 
stability if q = 0. 

Taking iaRe s (l7) +i(3(19) we get 
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iaResU(iaResu + i/3w) — iu(iaResu + i/3w) + U(iaResii x + i/3w x ) + 

+iaResU y v — f3 2 p + - — a 2 p = (-— — k 2 )(iaResu + i(3w). (28) 



Using the continuity equation (16), (28) becomes 

iaResU(u x + v y ) — iu>(u x + v y ) + U^-(u x + v y ) + 
-iaRe 5 UyV + k 2 p - = (— - k 2 )(u x + v y ). (29) 

Operating across (29) with d/dy, using the assumptions in boundary layer and 
the equation (18) we get the modified Orr-Sommerfeld equation for the boundary 
layer 

d 2 „ d 2 

(iaUReg — ica)(—z — k 2 )v — iaResU v = (— — — k 2 ) 2 v. (30) 
dy l dy l 

Considering normalization with viscous effects we would have 

(iaURes - - k 2 )v - iaRe s U v = (— - k 2 ) 2 v. (31) 

iL&g ay ay 

5 Stability Analysis 



We put equations (30)-(31) respectively in the form 

{Re s U{— - k 2 ) - Re s U" + - k 2 ))v = c{— - k 2 )v (32) 

and 



(Rep(-^ - k 2 ) - Rep" + ia-'Re 2 ^ - k 2 ))v = - k 2 )v. (33) 

In order to investigate the application of the Squire theorem, consider first the 
normalization with inertial effects (1321 . Often, however, we are only interested in 
the instability that appears first as the control parameter Re$ is increased. In this 
case, Squire's theorem tells us that we need only consider 2D disturbances. 

Consider a base state U(y). Imagine a growing 3-D disturbance to this base 
state at Reynolds number Re$3D, with wavenumbers a S o, Psd, and k\ D = a\ D + 
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(3 2 D . This corresponds to a solution c, v with c$ > 0( of the modified Orr- 
Sommerfeld equation 

d 2 d 2 
(ia 3D URe S 3D - ica 3D )(— - kl D )v - ia 3D Re S3 DU" v = {-^ - k\ D ) 2 v. (34) 

Now consider a 2D disturbance at a Reynolds number Re§2D- This has (3 2 d = 0, 
k2D = «2D and must satisfy the 2D modified Orr-Sommerfeld equation 

d 2 „ d 2 

{ia 2 DURes2D-ica2 D )(—-al D )v-ia 2D Res2DU''v = {-^-a 2 2D ) 2 v- (35) 

For values a 3 oRes 3 D = &2dRz&2D, «3D = ol 2 d and k 3D = a 2 D, this 2D modified 
Orr-Sommerfeld equation has the form 

d 2 „ d 2 

{ia 3D URes3D-ica 3D ){—-al D )v-ik 3D Re S 3 D U"v = (^~^3d) 2 ^> ( 36 ) 



which is exactly the same as [34J It must therefore have the same growing solution 
c, v with q > 0. 

Therefore, corresponding to the growing 3-D disturbance at Re$ 3 o with a 3 o, 
/3 3 d, and k 2 D = a 2 D + /3 2 D , there exists a growing 2D disturbance at Re^D with 
k 3 D = &2D and a.2D = &3D- These conditions lead to k 3 o = ol 3 d and we get 
that the disturbances are two-dimensional. Finally, we get Re^D = R&S3D, not 
R&&2D < R&53D- And so, we can not applied the Squire's theorem in the boundary 
layer of the flow. 

Using the same assumptions in (1331) . we get first Re 2 2D = ^^Re 2 3D with 
k 2 D = k 3D i.e. a 2D = k 3D and secondly Re 52 D = ^Re 2 S3D with k 2D = k 3D 
i.e. a 2 D = k 3D . So we have first Re^D > R^S3D and secondly Re^u < R&53D 
because k 3D > a 3D . We see therefore that we must take Re^D = R&63D and so 
k3D = &3D- We find the same result that the disturbances are two-dimensional 
and then Squire's theorem can't be applied. 

Finally, we still consider our three-dimensional disturbances so without using 
the theorem of Squire. Thus we write a = kcosO with 9 = (k x , k). This will 
allow us to deduce numerically whether the application of Squire's theorem in the 
boundary layer. Indeed 9 = corresponds to a two dimensional disturbance i.e 
k = a. 

We employ Matlab (Windows Version) in all our numerical computations to 
find eigenvalues. A Poiseuille flow with the basic profile 

U(y) = l-y 2 (37) 

is considered. 
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The eigenvalue problems (32)-(33) are solved numerically with the suitable 
boundary conditions. The solutions are found in a layer bounded at z — ±1 with 
U(±l) = 0. The results of calculations are presented in the following figures. In 
each group of figures, the first one is figure a), the second is figure b) and the third 
is figure c). First we present the figures relatives to the eigenvalue problem (32). 

For a fixed k = 1, we get figure 1 of q vs R for sequential values of 9 in which 
a) shows the entire graph, b) and c) are the magnified versions of a). 

For a fixed k = 2, we get figure 2 of q vs R for sequential values of 9 in which 
a) shows the entire graph, b) and c) are the magnified versions of a). 

For a fixed A; = 4, we get figure 3 of q vs i? for sequential values of # in which 
a) shows the entire graph, b) and c) are the magnified versions of a). 

Secondly we present the figures relatives to the eigenvalue problem (33). 

For a fixed k — 1, we get figure 4 of q vs R for sequential values of 9 in which 
a) shows the entire graph, b) and c) are the magnified versions of a). 

For a fixed k = 2, we get figure 5 of q vs R for sequential values of 9 in which 
a) shows the entire graph, b) and c) are the magnified versions of a). 

Through the figures 1 and 4, it is easy to see that if we take a curve with 9^0 
in the instability area (q > 0), we don't have necessary in the two normalization 
Re$(9 = 0) < Res(9 ^ 0). We also see through the figures that the normalization 
of time and pressure by inertial effects or viscous effects lead to the same order 
of stability /instability in the boundary layer. This confirms that in the boundary 
layer, viscous forces are on the same magnitude as the inertial forces i.e. the 
local Reynolds number is on unity order. We see also through the figures that 
at low Reynolds number the flow is stable but if the Reynolds number increase, 
unstability appears. The increase of the wave number induce also the stability of 
the flow. By figures 6 and 7 we find the first value of Res and a for which the first 
transition from stability to instability occurs i.e. q = in the case 9 = 0. We find 
the values (Res, a) = (5772, 1.02) which corresponds exactly to the one we know 
as critical value of the neutral stability of Poiseuille flow. Note that this value is 
the same in figure 6 — a) which corresponds to normalization by inertial effects 
and in figure 6 — b) which corresponds to normalization by viscous effects. We 
noticed also that for all overs values of k in the case 9 = correspond the same 
values of Res at q = whatever the normalization. We therefore conclude that 
in the boundary layer with a 2D-disturbance, we have the same neutral stability 
curve whatever the normalization. 
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Figure 2: (a) Growth rate (7, vs. Reynolds number Re for k = 2; (b) zoom of (a) 
to small values of k; (c) zoom of (b) to small values of Re. 
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Figure 4: (a) Growth rate Cj vs. Reynolds number Re for fc = 1; (b) zoom of (a) 
to small values of k; (c) zoom of (b) to small values of Re. 
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Figure 5: (a) Growth rate C, vs. Reynolds number Re for k 
to small values of k; (c) zoom of (b) to small values of Re. 



2; (b) zoom of (a) 



6 Conclusion 



In this paper, we have investigated the stability of boundary layer in Poiseuille 
flow. We have shown that the instability of the perturbed flow is governed by a re- 
markably equation named modified Orr-Sommerfeld equation. Contrary to what 
one might think, we find that Squire's theorem is not applicable for the boundary 
layer. We find also that normalization by inertial or viscous effects leads to the 
same order of stability or instability. For the 2D disturbances flow (9 = 0), we 
found the same critical Reynolds number for our two normalizations. This value 
coincides with the one we know for neutral stability of the known Orr-Sommerfeld 
equation. We noticed also that for all overs values of k in the case 9 = corre- 
spond the same values of Res at q = whatever the normalization. We therefore 
conclude that in the boundary layer with a 2D-disturbance, we have the same 
neutral stability curve whatever the normalization. 
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Figure 6: (a) Growth rate Cj vs. Reynolds number Re for k = 1.02 in the case 
of inertial effects normalization showing the critical Reynolds number for 2D- 
disturbances; (b) Growth rate Cj vs. Reynolds number Re for k = 1.02 in the 
case of viscous effects normalization showing the critical Reynolds number for 
2£>-disturbances. 



14 



References 

[1] J. G. Oldroyd. On the formulation of rheological equations of state. Proc. 
Roy. Soc. London Ser. A. 200 (1950), 523-541. MR 11,703a. 

[2] T. C. Ho and M. M. Denn. Stability of plane Poiseuille flow of a highly 
elastic liquid. J. Non- Newtonian Fluid Mech. 3 (1977), 179-195. Zbl 
414.76008. 

[3] Bertolotti, F. R, Herbert, T. and Spalart.R R Linear and nonlinear stability 
of the Blasius boundary layer.]. Fluid Mech, 1992, 242, 441-474. 

[4] K. C. Lee and B. A. Finlayson. Stability of plane Poiseuille and Couetteflow 
of a Maxwell fluid. J. Non-Newtonain Fluid Mech. 21 (1986), 65-78. Zbl 
587.76059. 

[5] R. C. Lock. The stability of the flow of an electrically conducting fluid be- 
tween parallel planes under a transverse magnetic field. Proc. Roy. Soc. 
London Ser. A. 233 (1955), 105-125. MR 18,700a. 

[6] Abu-Ghannam, B. and Shaw, R. Natural transition of boundary layers - 
the effects of turbulence, pressure gradient, and flow history. J. Mech. Eng. 
Sci.1980- 22 (5), 213-228. 

[7] Andersson, P., Berggren, M. and Henningson, D. Optimal disturbances and 
bypass transition in boundary layers. Phys.Fluids 11 (1),1999, 134-150. 

[8] Andersson, P., Henningson, D. and Hanifi, A. On a stabilization procedure 
for the parabolic stability equations. J. Eng. Math. 1998, 33, 311-332. 

[9] Landau, L. & Lifchitz, E. M. Fluid Mechanics. Butterworth Heinemann, 
Oxford, 1997. 

[10] Ray, Samad, and Chaudhury: Department of Mathematics, University of 
Burdwan, Burdwan-713104, West Bengal, India. 



15 



This figure "art.png" is available in "png" format from: 



http://arxiv.org/ps/1303.0535vl 



This figure "fig7.png" is available in "png" format from: 



http://arxiv.org/ps/1303.0535vl 



